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I investigate the propagator of the Wigner function for a dissipative chaotic quantum map. I show that a small 

amount of dissipation reduces the propagator of sufficiently smooth Wigner functions to its classical counterpart, the 

Frobenius-Perron operator, if ?i — *■ 0. Several consequences arise: The Wigner transform of the invariant density 

matrix is a smeared out version of the classical strange attractor; time dependent expectation values and correlation 

functions of observables can be evaluated via hybrid quantum-classical formulae in which the quantum character 

enters only via the initial Wigner function. If a classical phase-space distribution is chosen for the latter or if the 

map is iterated sufficiently many times the formulae become entirely classical, and powerful classical trace formulae 

0\ ' apply. 
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I. INTRODUCTION 



This is the third of three articles in which I describe semiclassical methods for dissipative chaotic quantum maps. 
The first two articles mM were devoted to spectral properties of such maps. It turned out that in the presence of a 
small amount of dissipation important information about the spectrum of the propagator can be obtained directly from 
the Frobenius-Perron propagator of the phase space density for the corresponding classical map. Spectral properties 
usually determine the behavior of most physical observables. It is therefore natural to ask what can be said about 
Qs ' expectation values of observables and their correlation functions, which are the most interesting quantities from an 
f~^ ' experimental point of view. 

(^ : In this article I lay a semiclassical framework that allows to calculate observables, correlation functions, and the 

^^ ' invariant state of certain dissipative quantum maps. As a central result it turns out that the propagator of the Wigner 

function is to first order asymptotic expansion in h identical to the Frobenius-Perron propagator of the phase space 

density for the corresponding classical map. Many important consequences follow. In particular, time dependent 

expectation values and correlation functions are given by quantum classical hybrid formulae, in which the quantum 

Ph ' character enters (to lowest order in h) only via the initially prepared Wigner function. If this function is a classical 

>^ ' phase space density, or after the map is iterated many times, the time dependent expectation values and correlation 

functions are given by entirely classical formulae. It follows that highly developed and precise classical periodic orbit 

Q ' theories can be applied. 

Ci There is currently very strong interest in understanding the interplay between chaos, dissipation and decoherence 

1-^ • in quantum mechanics p-[l5[. Dissipative quantum maps have been the object of choice in this area for a long time. 
The very same model as is studied here, namely a dissipative kicked top was investigated earlier numerically by 
Peplowski et ai, and I am going to derive some of their results analytically. Particularly noteworthy is the study by 
Graham and Tel on the quantization of the Henon map y . They also examined the time evolution of the Wigner 
function and found by different means very similar results for the propagator as will be derived in the present work, 
C^ ' namely a connection to the propagator of the entire dissipative classical map. Very similar results were also derived 
in the theory of superradiance, where the quantum character of the problem enters in an initial distribution of points 
in phase space, each then giving rise to a classical trajectory. In this way the macrosopically amplified quantum 
fluctuations of the delay time of the superradiant pulses was predicted |l^ . 

The paper is structured as follows. After introducing in the next section the type of quantum maps that will be 
dealt with in this paper, I will derive in section fl| a semiclassical approximation of the propagator of the Wigner func- 
tion. Section |^ exploits the consequences for the invariant state, expectation values of observables, and correlation 
functions. The results are summarized in section M. 

II. DISSIPATIVE QUANTUM MAPS 

A dissipative quantum map P is a map of a density matrix p from a time t to a time t + T: 

p{t + T) = Pp{t). (2.1) 

The density matrix should be thought of as a reduced density matrix, resulting from a larger one, with the degrees 
of freedom of the environment which causes the dissipation traced out |17[ . The type of maps that I consider in the 



following have been described in much detail in the earlier papers KM ■ I will therefore restrict myself to introducing 
them here only very briefly, referring the interested reader to the above references. 

The maps are particularly simple in the sense that the dissipation is well separated from a remaining purely unitary 
evolution where the latter by itself is capable of chaos. The unitary part is described by a unitary Floquet matrix 
F, and the dissipation by a propagator D. The prime example will be a dissipative kicked top, which is an angular 
momentum J whose orientation alters due to a continuous time evolution and periodic kicks. Between two unitary 
evolutions the angular momentum experiences damping. The dynamical variables of the top [ |l8| , pi)[ are the three 
components Jx.y,z of the angular momentum. Consider a unitary evolution generated by 

p^Q-iijJ^Q-i^Jy ^ (2.2) 

This Floquet matrix first rotates J by an angle /3 about the y-axis and then subjects it to a torsion, i.e. a rotation 
with a rotation angle proportional to J^ about the z-axis; the torsion strength is given by the parameter k. The 
maps considered (including the dissipative parts) have the important property that the square of J is conserved, 
J^ = j{j + 1) = const, with j a fix ed p ositive integer or half integer. They also have a well defined classical limit 



formally attained by j -^ cxd. In (2.2) I have introduced J = j + 1/2 which simplifies most of the semiclassical 



formulae. It can be shown that the Bloch sphere limj^oo J^/J^ = 1 is the classical phase space [|19|; ^i — coa6 plays 
the role of classical momentum and 4> the role of the conjugate coordinate. The angle 6 denotes the polar angle of 
J reckoned against the Jz axis, the azimuthal angle. The dynamics generated by F alone has been extensively 
studied in the literature ||l^,p9| and is known to become strongly chaotic for sufficiently large values of k and /3. The 
parameter values fc = or /3 = lead to integrable motion. 

As damping mechanism I consider a process which is given in continuous time r by the Markovian master equation 

^Pir) = ^([J-, P(r) J+] + [J-p{r), J+]) ^ Ap(r) . (2.3) 



The linear operator A is hereby defined as generator of the dissipative motion. Equation ( p. 3D is well-known to 
describe certain superradiance experiments, where a large number of two-level atoms in a cavity of bad quality 
radiate collectively [p0[[2l]] . The angular momentum operator J is the n th e Bloch vector of the collective excitation 



and the J+, J_ are raising and lowering operators, J± = Jx ^ iJy Eq.(2.3) is formally solved by p(r) = exp(AT)/7(0) 



for any initial density matrix p(0), and this defines the dissipative propagator 

D{t) == exp(AT) . (2.4) 



Expficit forms of D can be found in |20|j2^,|23|l . 

Damping manifests itself in a reduction of the J^ component which in the quantum optics application measures the 
energy stored in the two-level atoms. In the classical limit, formally attained by very large values j, the Bloch vector 
creeps towards the south pole 6* = tt of the Bloch sphere. The corresponding classical maps for rotation, torsion and 
dissipation can be found in the appendix of |^. The time r is measured in units of the classical time scale. In the 
following r will be set equal to the time between two unitary steps, and therefore measures the dissipation strength. 
The total time will be measured in units oi T, t = NT, and I will only keep track of the discrete time N. 



Given F by (2.2) and D by (2.4), the total map reads 

p{N+l) = D{Fp{N)F'') = Pp{N). (2.5) 

I have suppressed the dependence on the system parameters fc, (3 and r. Due to the dissipation the total propagator 
P is non-unitary. It has always one eigenvalue equal to one which corresponds to an invariant density matrix. Its 
existence follows solely from probability conservation ||l^. All other eigenvalues have absolute values smaller than 
one, reflecting the dissipative nature of the map. Correspondingly, the classical map leads to a shrinking phase space 
volume and, in the case of chaos, typically to strange attractors. 

In the following a semiclassical approximation for P will be of importance which has been derived in M . I write it 
in the Jz basis {Jz\j,m) = m\j,m)) and use indices m and k related to mi, m2 of pmim-iiN) = {j,mi\p[N)\j,m2) = 



Pm{k, N) by ?7i = (jni + TO2)/2 and k = (mi — m2)/2. In such a representation (2.5) reads 



Pm{k,N+l)=Y,Y. Pmk-m'k'Pm'{k', N) (2.6) 

7n' k' 

I convert the discrete sums into integrals by Poisson summation. For large J it is convenient to go over to rescaled 
coordinates, p — m/ J and 77 — k/ J; the density matrix p,„(fc, N) will then be denoted by p{p, rj, N). Let us think of 



it as a continuous function of the continuous variables /x and 77. Denoting the summation variables from the Poisson 
summation by si and S2, the new density matrix p(/z, ry, iV + 1) in continuous arguments is obtained from the old one, 
p(/x, 77, N) according to 

/oo 
d^l'dr^' Yl e'2'^'^(^^(^'+''')+^^(^'-"'))p(Ai,77;/^',7y')p(M',^',A^)- (2.7) 

si .S2 — — 0C 

Note that the Poisson summation was done in the original quantum numbers m'l , m2 in order to avoid problems 
arising from the fact that m' and k' can be half integer. The factor two in front of the integral arises from the back 
transformation to m' , k' and thus jjl' ,ri' . 

Central object of interest for the subsequent study is the total propagator P. A semiclassical approximation has 
been derived in great detail in ||l|. Here I just summarize the main features that will be of importance in the following 
and refer the reader interested in the details of the derivation to p| . 

The semiclassical form of P is very reminiscent of a van Vleck propagator p4| . It is given by a double sum over 
all classical paths ai and CT2 joining the initial and final coordinates (/i', cj)') and (/i, (p), respectively. The double sum 
arises because we propagate a density matrix and not a wave function. Along every path a complex "action" V' is 
accumulated. It contains a real part R from the dissipative part of the motion and two imaginary components iS^-^ 
and —\Sa2 from the Floquet matrices F and F\ 

ipifi, ri\ /i', 7]') = R{ii, V; 77) + 1(5*^1 (v + 77, ^' + 77') - Sa2 [T^ -r],fi' - 77') + 2TrW) . (2.8) 

The dependence of ip on its arguments arises from its implicit dependence on the intermediate coordinate 77 = 
F(/i, 77; /i', 77') which is the solution of 

dj7ip = 0. (2.9) 

The integer I stems from a Poisson summation that is still present in the propagator displayed below. The propagator 
also contains three pre-exponential factors, a function B from the dissipation, a C^-^ from F and another C*^ from 
Ft. They are basically second mixed derivatives of the actions, 

C,{^,^) = ^^J\d,d^S,\, (a = ai,a2), (2.10) 

Both derivatives in the last equation have to be taken at constant 77, and the second one additionally at constant 
energy E of a fictitious Hamiltonian system that underlies the dissipative dynamics |^^. These distinctions will be 
of no further importance in the remainder, though, as the saddle point approximations that will come up soon, pick 
classical trajectories at 77 = 0. For these both derivatives under the square root are the same and both combine to 
the Jacobian connected with the classical trajectory [|l|. 
In its full-fledged form the propagator reads 

oc 

P(/i,77;/i',77')= Y^ ^ ^B(M,J7;77)ai(l7 + 77,Ai' + ?7')C(^-7?,Ai'-7?') 

i — — CXD (y\ ,(72 77 

exp(j7/>(/^,77;/^',77')) . (2.12) 

The explicit form of the actions S and i? is not relevant here. More important are their generating properties for 
the classical trajectories, 

df,Sa{l^,fJ.) ^<Pl{'^,fJ.), duSa{v,^i)^-(j)i{v,^i), (2.13) 

where 0^ and 0^ are the initial and final coordinates of the classical trajectory a corresponding to the unitary part 
of the map {a — ci, (72). For R we have 

^^,R{^l, I/; 0) = ^ /i = ^ld{v) ^ d,,R{fi, iy;0) = 0, (2.14) 

where /i — fJ-di^) denotes the (unique) classical trajectory corresponding to the dissipative part of the motion. Fur- 
thermore, i? as a function of 77 has a single maximum at 77 = 0, i.e. 

5,,i?(Ai, I/; 77) = ^ 77 = . (2.15) 



III. WIGNER FUNCTION AND WIGNER PROPAGATOR 

In order to unravel classical properties of the quantum map, it is natural to go to a phase space representation of 
the density matrix. It turns out that the Wigner function is very well suited for this purpose. In fact, the Wigncr 
function has been used many times in order to study the transition from quantum to classical mechanics 0-0,0,121 . 
I will show in this section that the propagator of the Wigner function is — for sufficiently smooth Wigner functions 
— nothing but the classical Frobenius-Perron propagator of the phase space density. 

Usually the Wigner transform is defined as a Fourier transform with respect to the skewness of the density matrix 
in coordinate representation [fiSpS], 

1 1"°° T T 

(P,9) = ^y ^rf^e'^-/''(g--|p|g+-). (3.1) 



Pw 



In our problem we have p in the momentum basis /i. Inserting resolutions of the identity operator in the momentum 
basis in the above definition of pw (p, q) we obtain 

Pw{p,q) = Tr^ dee''«/''(p+||pb-i). (3.2) 

Note the change of sign in the skewness. We have p = p,, q = (j>, and 1/ J replaces Ti. An additional factor J arises 
because the original quantum numbers m and k are rescaled to p and rj as explained above. I therefore define the 
Wigner transform of p{p, rj, N) as 



p^(p,^,N)^—J drje''^''''p{p,j,N). (3.3) 

It has the right normalization in the sense that 

/ dpd(j,pw{p,(l),N) = J / dpp{p,0,N) ~ ^p™™(Ar) = 1. (3.4) 

ra 

Note that the corrections from passing from the integral to the discrete sum of the diagonal matrix elements are of 
order 1/ J and become negligible in the limit of large J, as long as p{p, 0, N) does not fluctuate on a scale 1/ J, i.e. as 
long as the probability profile has a classical meaning. 
The inverse transformation reads 

p{p, Tj,N)^jJ d</)e-'2^"*pv^(M, ^, N) . (3.5) 

Wigner functions on SU{2) have been introduced before in the literature |pq,P2l ^^^ angular momentum coherent 
states and appropriate transformations of Q- or P-functions. These definitions avoid problems at the poles of the 



Bloch sphere that can arise in the present approach. On the other hand, the definition (3.3) is much simpler from a 
technical point of view and sufficient for our purposes. 

We are now in the position to calculate the Wigner function after one a ppli cation of the map from the original one. 



To this end we insert the propagated density matrix p[p, rj, N + 1) from (2.7) in 



t2 /-oo 

Pm/(m>,^ + 1) = — / drje^'J^^p{p,rj,N + l), (3.6) 

7^ J-oo 

and then express the original density matrix p{p' , rj' , N) in terms of its Wigner transform, 

pwifJ-, (t),N +l) = / drjdp' dri' dcj)' ^ exp (2iJ(7r((si + S2)p' + (si - sa)??') - Tl'<i>' + V'f))) 

Sl,S2 

P{p,mp',rj')pw{p\4>\N). (3.7) 

With the semiclassical expression ( ^.12| ) for the propagator, we arrive at 



2 j3 r 

pw{P;(f),N +1)^ — I dr]dfi'dr]'d(j)' ^ B{fi,V;r])Ca^{V + ri, fi' + ri')C*^{V - r], f/ - r]') 

exp(JG(/x,7y;^',7/)) pw{fi\^\N), (3.8) 

where the "action" G is given by 

G{p, r,; p', r,') = ^{p, r,; p\ r^') + 2t^cJ) - 2r,V' + 2tt{s^{p' + r/) + s^^/ - v'))) ■ (3.9) 

The form of the integrands and the fact that P is already approximated semiclassically, i.e. correct only to order 1/ J, 
suggests to integrate by saddle point approximation (SPA). To do so, we must assume that the initial Wigner function 
pwip-' , <t') is sufficiently smooth, i.e. has no structure on a scale 1/J. This is at the same time a necessary condition 
if we want to attribute a classical meaning to pw- 
The saddle point equations read 



a„G = d,,R + chjRdrjV + i ((-0^^ + q^i^ + 2ttI) dnV 
df^^G = djjRdf^.V + i ((-0^^ + q)L + 2t:1) d,,V 

d^,G = 2iri' = Q. 






""1 



2ttI^ dri'V 



0^^ -0^^+20)^0, 
(j)l^_ + 27r(si + S2)) - 
</.;^ + 27r(si - S2) - 2^') 



(3.10) 
(3.11) 
(3.12) 
(3.13) 



For brevit y I ha ve suppressed the arguments of R and (f/^jCJjl. They are the same as in ( p.8[ ) for R and So-, cr = cti, (72. 
Equation ( ^.13| ) immediately gives 77' = 0. To solve the rest of the equations, let us first assume that d^jR = 0. I 
will show below that this is the only possible choice. Then we have from the generat ing p roperty (2.14) that p, is 
connected to V via the classical dissipative trajectory, p = pd{i^) and the real parts of (3.11) and ( [3.1^ )~give already 
zero. I will assume that all relevant solutions to the saddle point equation be real, as is expected from the physical 
origin of the variables as real valued quantum numbers. Real- and imaginary parts of all saddle point equations must 
then separately equal zero, so that we confront eight instead of four equations. The assumption djjR — solves two 
of the m at the same time. The real part from (3.10) gives additionally 9,,i? — and thus according to the property 
(2.15) 77 = 0. Only the propagation of probabilities, i.e. the diagonal elements of the density matrix contributes in 
the saddl e po int approximation; the same was true for the calculation of trP^ |I| . 

From (2^) follows i (— ^^i + ^^2 + 27r/) = 0, i.e. —4>l^{T',p) + (j)l^{u,p) + 2ttI — 0. Thus, the final canonical 
coordinates of the two trajectories ui and (72 raust agree up to integer multiples of 27r. Since also initial and final 
momenta are the same {p and V, respectively), the two trajectories must be identical, i.e. cti = 0-2 = <7. Counting all 
angles modulo 27r we also have / = 0. 

The imaginary part of (3.10) leads to (p^ = (p, the imaginary part of (3.11) to ,si + S2 = 0, and the imaginary part 
of ( 3.12| ) to (pl^ = (j)' + 2t:s2- These equations describe precisely the classical trajectories for the unitary part of the 
motion from an initial phase space point [p' , (/)') to a final one {y,(t>), again counting the angle modulo 27r. Together 
with p ~ Pd{p) the saddle point equations thus give the classical trajectory from (/i',(/)') to {p,(t>). Note that this 
trajectory is unique if it exists, since classical trajectories are uniquely defined by their starting point in phase space. 

For evaluating the SPA we further need the determinant of the matrix G*^^-' of second derivatives of G. It is 
straightforward to verify that its absolute value is given by 



|detG(2)| = 16|5i7C(^:M')l' 



(3.14) 



The overall phase arising from the SPA equals zero, as can be seen by the same techniques that were used in [n|. 
Calling the classical map / : (/i', 0') -^ {p,<P) ~ f{p'i4'')i we get the saddle point approximation 



pw{p,(f>,N + 1) 




B{pdiv),V;0)\CiV,p')\^pwir\l^,^),N) 



Pw{r'{p,^),N). 



(3.15) 
(3.16) 



The prefactor in the last equation is nothing but the inverse of the Jacobian of the classical trajectory which arises 
solely from the dissipative step since the unitary one has Jacobian unity. So with the abbreviation y = {p,4'), 
X — (p' , (f)') of final and initial phase space coordinates we have 



pw{y,N + l) 



Pw{f-Hy),N) 



a/ 

dx 



dx5{y - f{x))pw(x,N) = I dxPw{y,x)pw{x,N) 



(3.17) 



^=/~Hy) 



which identifies the propagator of the Wigner function as the classical Frobenius-Perron propagator of the phase 
space density, 

Pw{y,x) ^ Pciiy,x) ^ S{y - fix)) . (3.18) 

Note once more that this conclusion holds only if the test function pwi^) on which Pw acts is sufficiently smooth, 
namely does not contain any structure on a scale 1/J or smaller. For classical phase space densities this is often not 
the case. Indeed, continued application of a chaotic map leads to ever finer phase space structure, so that after the 
Ehrenfest time of order A~ ^ In J (where A is the largest Lyapunov exponent) scales are reached that are comparable 



with 1/J. From equation (3.18) it follows immediately that also the propagators of the iterated maps are identical, 

W 



P^ = Pjy, but the validity is restricted to discrete times N much smaller than the Ehrenfest time. 



Let me finally show that there is no alternative to the assumption djjR = about the solution of the saddle point 



equation if only real so lutio ns are permitted. To see this suppose that djjR ^ 0. Then we have from (3.11) that 
dv/dfj,' = 0, and from ( [3.12[ ) dv/drj' = 0, such that 77 is a function of ^ and 77 alone. The imaginary part of ( 3.11 ) 
gives 0^^ — 0^2 + 27r(si + S2) = 0. If we differentiate with respect to /i' and with respect to 77' and remember that 
ry' = follows directly from (|3.13|) , we are immediately lead to d^i(j/^_^ {y + 77, /i') = d^icfjl^^ {V — -q, /i') — 0. Thus, all 
trajectories with given initial 0^ end at the same final momentum u -\-rj (for a — ai) or V ~ rj (for a — 02). From the 



imaginary part of (3.10) follows in the same fashion dp,i(^{,^ {p + rj, /i') = d^'(j)l^ {V — rj, //') = 0. So the final canonical 
coordinate does not depend on the initial momentum either. In other words, all trajectories with the same initial 0^^ 
(respectively (/)^ ) but arbitrary initial /i' end at the same final phase space point. But this is in contradiction with 
the fact that a final phase space point uniquely defines a trajectory. Therefore, the initial assumption djjR =/= must 
be wrong, q.e.d.. 

IV. CONSEQUENCES 



Equation ( 3.18| ) is a key equation from which many consequences follow in a straightforward way. Let me first show 



that previous results about spectral properties are readily recovered. 

A. Spectral properties 

In ||l|,g| we have shown with great effort that 



trP'^ = trP^i (4.1) 



for all integer iV, if J — > 00. The same result can now be obtained much more easily. In view of (3.18) all that 

W 



remains to do is to show that trP^ = tiPyy. To see this, let us extract the general relation between any Pw and the 



corresponding P from (3.7) and the definition of Pw in ( ^.17 ). Comparing the two equations we are lead to 

Pw(A^>;/i',0') ^— [ dT^dri'Y e2'^(''>'-''-^)+'2-4(^i+^^)'^+(^i-^^)''')p(;,,ry;/i',7y'). (4.2) 



Sl,S2 



The equation holds for any propagator P of the density matrix, therefore also for the propagator P of the A^th 

^N 



iteration of the original map. It is then one line of calculation to show that the trace of P^ 



trP^ =y d/z#P^ (/.,</.; /i,0), (4.3) 

is given by 

trP,!^ = 2J2 fdT]d^l J2 e''"''«^^+"^)^+("i-"^)'')P^(M,?7;M,?7) =trP^, (4.4) 



Sl,S2 



where in the last equation I have gone back to discrete summation, undoing the Poisson summation. Thus, we have, 
up to 0(1/ J) corrections, trP^ = trP^ — trPj^. The easiness of the derivation comes along with a more severe 
restriction on the validity of the proof, though. As mentioned, we need A'^ ^ A~^ln J, whereas in the much more 
involved calculation avoiding the Wigner function H,0] the discrete time N had only to be much smaller than the 
Heisenberg time of the non-dissipative kicked top, TV <C J. 

6 



B. The invariant state 

If the classical and the Wigner propagator are the same up to order 1/ J corrections, so are their eigenstates. An 
invariant state of P is defined as an eigenstate with eigenvalue one, Pp(oo) — p{oo), and correspondingly PwPw{oo) — 
Pvf(oo), PciPci{oo) = Pciioo). Classically, there can be many invariant states. For example a distribution Pci{x) — 
d{x — Xfp) is an invariant distribution if Xfp is a fixed point in the neighborhood of which the map preserves area. If 
several fixed points exist one can linearly combine the delta functions on them to obtain new invariant distributions. 
In the present context we are interested, however, in invariant distributions that are not only stationary, but can also 
be obtained from a general initial state by iterating the map infinitely many times. I indicate this with the arguments 
"infinity". If the system is ergodic, the final invariant distribution is unique. For volume preserving chaotic maps, it 
is a homogeneous distribution in phase space regions selected by the remaining integrals of motion. For dissipative 
chaotic maps one encounters typically a strange attractor in phase space psf. 



From ( 3.18 ) we conclude that up to 0(1/ J) corrections 

pwi'x) ^ Pcii'x) . (4.5) 

The corrections have to be understood as a smearing out on a scale 1/ J. Indeed, suppose we start from a smooth 
initial Wigner function, and then iterate it many times with Pw, it evolves according to (|3.18) up to the Ehrenfest 



time as a classical phase space density. After the Ehrenfest time the classical dynamics continues to produce ever finer 
structures in the phase space den sity, whereas Heisenberg's uncertainty principle prohibits structures of pw smaller 



than 1/ J. As pointed out before, ( p_J-8[ ) therefore ceases to be valid, and pw is left at the stage where it is the smeared 
out classical strange attractor. Figjl] shows that indeed the quantum strange attractor is a smeared out classical one. 
The Wigner function was obtained by direct diagonalization of the propagator and subsequent Wigner transformation 
of the eigenstate with eigenvalue unity, the classical picture by iterating classical trajectories many times and making 
a histogram in phase space. Similar observations were made earlier numerically P Jl^ , |29| , |30t . 

C. Expectation Values 

Suppose that a system is prepared at time i = by specifying the density matrix p(0), or equivalently the initial 
Wigner function pw{x,0). We let the system evolve for a discrete time N and then measure any observable A of 
interest. The expectation value of the observable is given by 

{A{N)) = tr(ip(7V)) = I dxAw{x)pw{x,N) , (4.6) 

where Aw{x) is the Weyl symbol associated with the operator A. The definition of A is analogous to the definition 
of Pw |1^- So Aw{x) is also a phase space function. To lowest order in 1/ J it equals the classical observable A{x) 



that corresponds to A, if the classical observable exists. Using (3.18) we immediately obtain 



{A{N)) = / dxA{x)P^iPw{x,Q), (4.7) 

up to corrections of order 1/ J. Thus, quantum mechanical expectation values can be obtained from the knowledge 
of the classical propagato r an d the classical observable for any initial Wigner function that contains no structure on 



the scale 1/J. Equation (4.7) is a hybrid classical-quantum formula, since the initial Wigner function can be very 
non-classical, e.g. can contain regions where pw{xtG) < 0. 

But also the cases are interesting where already the initial Wigner function is a classical phase space density, 
Pw = Pch or where the time N is sent to infinity, such that the invariant state is reached. In both cases the quantum 
mechanical expectation value is given by a purely classical formula. In particular, expectation values in the invariant 
state pw{oo) are given by 

{A)^= / dxA{x)pci{x,(x) , (4.8) 



since up to order 1/J corrections PdPwioo) = PciPci{oo) = pci{oo). This allows us to use highly developed classical 
periodic orbit theories to evaluate {A)cx> p^]. These theories can be made very precise by the so called cycle expansion, 
which means that classical prime cycles p with similar actions are systematically combined to so-called pseudo orbits 



or cycles tt. To expose these theories would be beyond the scope of this article. Let rae just present the result for {A)oq, 
explain how one uses it and show numerically that the agreement with ab initio quantum mechanical calculations is 
indeed very good. 

Starting point for the practical use of the classical trace formulae is a list of prime cycles of the classical map, their 
stabilities, and topological lengths which has to be calculated numerically. Prime cycles are periodic orbits that can 
not be divided into smaller periodic orbits. The stabilities of a prime cycle are the (in the present context: two) 
stability eigenvalues, i.e. the eigenvalues of the Jacobian connected with the map from the starting point of the cycle 
to the last point before the cycle closes. Since phase space volume is not conserved for dissipative maps the product 
of the two eigenvalues usually does not equal unity, so we need to calculate always both of them. In the following, Ap 
will denote the product of all expanding eigenvalues (i.e. with absolute value larger than one) of a prime cycle p, and 
1, if both are contracting. The topological length is for maps just the length in discrete time, with the convention 
that fixed points of / have topological length Up — 1. Finally, we need the values Ap of the observable averaged along 
the prime cycles. Out of the prime cycles one has to construct all distinct non-repeating combinations {pi, . . . ,pk} 
with a given total topological length n^ = Up-^ + . . . + Up^ . The other characteristics of the prime cycles are combined 
to corresponding quantities for the pseudo cycles as well, At^ = Ap-^ + . . . + Ap^, , and Att = Ap^ • . . . • Ap^ ; and wc define 
the pseudo-cycle weight i^ = (— l)'^+^/|A7r|, where k denotes the number of prime cycles involved. With all this, the 
expectation value of the classical observable A in the invariant state, (A)oo = / A(x)pci{x, oo) dx is given by [ p2| 

(A)^ = ^^. (4.9) 

The advantage of the cycle expansion is that the periodic orbit sum is truncated in a clever way, such that different 
contributions that would lead to a highly fiuctuating behavior almost compensate. Fig.0 shows a comparison of 
exact quantum mechanical resul ts fo r the observables Jz/ J and Jy/ J in the invariant state, compared to results from 
classical periodic orbit theory ( |4.9| ), as well as straightforward classical evaluations. The latter were obtained by 
iterating many randomly chosen initial phase space points and averaging over the generated trajectories. Whereas 
Jy fluctuates only slightly about the value zero suggested by the symmetry of the problem, Jz/J decreases from 
at zero dissipation to —1 for strong dissipation as the strange attractor shrinks more and more towards towards 
the south pole of the Bloch sphere . The figure shows that even with rather short orbits (n^ < 4) the quantum 
mechanical result is produced very well for both observables. The agreement improves as expected with larger values 
of J and becomes almost perfect when comparing the classical simulation with ( [4.91 ). Instead of ordering the cycles 
by topological length, I also tried stability ordering 133], but did not observe further significant improvement. 

D. Correlation functions 

The discrete time correlation function K{N2, Ni) between two observables A and B with respect to an initial density 
matrix p(0) is defined as p3 

K{N2,Ni) = {B{N2)AiNi))Q = tr (BP^MP^V(O)) • (4.10) 

This function has typically a real and an imaginary part. The latter is connected to the Fourier transform with respect 
to frequency of a linear susceptibility (see e.g. u^)- Here I show how the real part of K{N2,Ni) can be calculated 
semiclassically. 



Starting point is the observation that in the above expression ( 4.10 ) AP^^p{0) ~ Ap{Ni) enters in the same way, 
as the initial density matrix p(0) enters in the calculation of {A{N)) (eq. |4.6| ). In fact, formally K{N2, Ni) is nothing 
but the expectation value of B with respect to the N2 times propagated "density matrix" p'{Ni) = Ap{Ni). Note 
that p'{Ni) is not really a density matrix, in general not even a hermitian operator. However, in the derivation of 
expectation values the special properties of a density matrix (besides hermiticity also positivity, and trace equal unity) 
did not enter at any point. The only thing that did matter was that the density matrix had to have a smooth Wigner 
transform. This I will suppose as well about the Wigner transform p^(iVi) of p'{Ni). Later on we will see that the 
assumption is justified if the initial density matrix p(0) has a smooth Wigner transform and the Weyl symbol Ayy 
a smoot h cl assical limit. So let us introduce a Wigner transform p'y^,{x,Ni) in complete analogy as for any density 
matrix (3.6), 



p'wi^^,<l>) = — / d77e2'^''^p'(^,,y,iVi) (4.11) 

n 



and then use ( [4.7D to express the correlation function K{N2, Ni) as 

K{N2, iVi) = [dxdy B,i{y)P^'p'^{x, iVi) 



Now I write p'{^i,ri,Ni) = {m + k\Ap(Ni)\m - k) in ( [4.1 1| ) as 

/oc 
dA ^ (TO + fc|yl|JA)(JA|p(iVi)|m-fc)e- 



i J27rnA 



(4.12) 



(4.13) 



n— — oc> 



where I have introduced a factor unity with Z — JX as summation variable and then changed the sum to an integral 
over I by Poisson summation. In terms of the corresponding Weyl symbol and Wigner function we have 



{m + k\A\l) 



1 

2^ 



{l\piN,)\m-k) = j 



l(t>ie-iip{-iJ{p + ri- X)(l)i)Aw{ >i) 

h exp (-« J(A - p + vi)cj)2) pw{ ;:; ,02, ^i) 



(4.14) 
(4.15) 



27r2 



E/ 7^ 7 11 11 A ,^J' + 1^ + ^ 1 \ A + fJ'-v 1 «r 
/ dXd-qd(j}id(t)2Aw{ ,(Pi)Pw\ 7. ,02,^1; 



n— — 00 



If we insert the last two equations into ( 4.13 ) and the resulting one into ( 4.11 ) we are lead to 

p'nr{p,(l),Ni) = 

exp(iJif(A,77,0i,(/)2)), 
with an exponent H given by 

H{X,ri,(j)i,(j)2) = 2?70- ifi + ri~ A)0i + (/i - 77 - A)(/)2 + 27rnA . 
Integration by SPA leads to the saddle point equations 

d^H = 20 - 01 - 02 = 

dxH = 01 - 02 + 2?™ = 
d^,H = -{p + T^-X) = 
9^2 H — p — rj — X = 0. 



(4.16) 
(4.17) 



(4.18) 
(4.19) 
(4.20) 
(4.21) 



The second one gives immediately 0i = 02 mod27r; and if we restrict 0i, 02 as before to a 2tt interval, we have n = 
and 01 = 02 = from (4. IS). The last two equations give ry = and X = p. The value of H at the saddle point is 
zero and one easily checks that the determinant of second derivatives gives 4. Putting all pieces of the SPA together 
we obtain 



p'yyix,Ni) = Aw{x)pwix.,Ni) . 



(4.22) 



This means that p\^{x,Ni) is smooth if Aw{x) and pw{x,Ni) are smooth. If we remember th at to lowest order in 
1/ J the Weyl symbols Aw and Bw are just the classical observables A and B, we obtain from (4.10) the final result 



KiN2,Ni 



dxB{f'^-{x))A{x)pw{x,Ni) 



(4.23) 



So semiclassically, the correlation function has the same structure as a classical correlation function with respect to 
a classical initial phase space density Pci{x), 



K,iiN2,m)^ / dxBif"-{x))Aix)pciix,Ni) 



(4.24) 



The only quantum mechanical ingredient is the Wigner function after A^i steps. We have the same kind of hybrid 
classical-quantum formula as for expectation values. And as for expectation values, in the limit of large A^i and with 
A^2 — Ni — N kept fixed, the quantum mechanical correlation function approaches its classical value, as pwix,Ni) 
tends to the smeared out classical invariant state Pci{x, 00). Nevertheless, as pointed out in the context of expectation 



values, pwix,Ni) can describe very non-classical states as for instance Schrodinger cat states pll 



Remar kable about ( 4.23| ) is also the fact that the expression is always real. We can trace this back to the realness 



of p'yy in ( 4.22 ). Since Ap{Ni) is not necessarily hermitian, there would be no need for p\^ to be real. However, note 
that Ap{Ni) would be hermitian, if A and p(A^i) commuted. Since they do commute classically, the commutator must 
be of order 1/J, and the imaginary part in K{N2,Ni) is therefore always at least of one order in 1/J smaller than 
the real part. 

If A^i ^ oo (with N — N2 — Ni fixed) or if p{x,0) is chosen as the invariant density matrix so that K{N + 
Ni,Ni)ni-*oc = K{N) = Kci{N) = Kci{N + Ni,Ni)]y^^ao, we can use classical periodic orbit theory to calculate the 
quantum mechanical correlation function ||35| . The use of the theory is completely analogous to the case of expectation 
values. In fact, the classi cal correlation function is nothing but the expectation value of B{N)A{0) in the invariant 



state pci{oo), so that in (4.9) we just insert for Ap the variable B{N)A{0) averaged along the prime cycle p. The 
practical evaluation of K{N) via the periodic orbit formula is, however, handicapped by the fact that for K{N) one 
should have at least prime cycles of the length N. Finding all of these for large iV is a difficult numerical problem, 
and hindered in our example additionally by the fact that we do not have a symbolic dynamics for the dissipative 
kicked top. Nevertheless, fig.(H) shows that at least the classical result and the real part of the quantum mechanical 
correlation function (./^(iV) J^(0)) agree rather well. 

V. CONCLUSIONS 

In this article I have shown that in the presence of a small amount of dissipation (r > 1/J, where J — > cxd in 
the classical limit), the propagator of Wigner functions that are smooth on a scale 1/J agrees up to order 1/J 
corrections with the classical Frobenius-Perron propagator of the phase space density. From this key result a number 
of important consequences followed. First, it allowed in a much simpler way than before to prove tvP^ = trP^ for 
fixed N as J —> 00. Second, it gave rise to compact semiclassical formulae for expectation values and correlation 
functions of observables. Basically, expectation values and correlation functions of observables can be evaluated with 
hybrid classical-quantum formulae, with the only remainder of quantum mechanics lying in the initial condition, 
i.e. the initial Wigner function has to be used to average over phase space instead of a classical phase space density. 
If one iterates the map sufficiently many times, the Wigner function approaches the smeared out classical strange 
attractor and from then on the difference between quantum mechanics and classical mechanics is no longer visible 
in expectation values or correlation functions. Powerful classical periodic orbit theories can be applied to address 
questions in quantum mechanics. Another way to obtain classical results is to start from the very beginning with a 
Wigner function that is a classical phase space density. 
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FIG. 1. Wigner function p„(oo) corresponding to the invariant density matrix {j — 40) and classical stationary probability 
distribution pci{oo) (strange attractor) for k — 4.0, /3 — 2.0, r = 0.5. The Wigner function is a "quantum strange attractor", a 
smeared out version of the classical strange attractor. The range of arguments is </!> = — vr . . . tt and /i = ~1 (in the background) 
... 1 (in the foreground). 
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FIG. 2. Three different ways of calculating expectation values in the invariant state as function of the dissipation: Quantum 
mechanically {j = 20, full lines), direct classical simulation (dashed lines), and classical periodic orbit formula with cycle 
expansion (circles: Jz/J, squares: Jy/J). The system parameters are fc = 8 and P — 2. 
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FIG. 3. The correlation function {J^{N)J:,{0))/J^ in the invariant state as function of A^ for fc = 8.0, /3 = 2.0, r = 1.0. The 
real part of the quantum mechanical function {j — 40, circles) agrees well with the classical result (dashed line). For clarity 
the classical result has been plotted for continuous N, even though it is only defined for integer A*'. 
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